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We show that the propagation of a TV-photon field in space and time can be described by a 
generalized Huygens-Fresnel integral. Using two examples, we then demonstrate how familiar Fourier 
optics techniques applied to a TV-photon wave function can be used to engineer the propagation of 
entanglement and to design the way the detection of one photon shapes the state of the others. 
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In quantum optics and quantum information science, 
one often deals with systems in which the total pho- 
ton number TV is known, the challenge being to prepare 
the required iV-photon state and engineer its evolution. 
When N is small, a wave-function formalism provides a 
more compact and intuitive physical description than the 
usual quantum field formalism. This has been advocated 
in the recent years [1, 2] and justifies the renewed interest 
in photon wave-mechanic. 

Single photons propagate in space and diffract on ob- 
stacles exactly as classical waves do. How a TV-photon 
wave function propagates is less obvious. In Ref. 1, the 
authors used a two-photon Maxwell-Dirac equation to 
study the disentanglement of a photon pair. In this 
Letter, we would like to lay the grounds of a propaga- 
tion theory of TV-photon wave-packets. Instead of using 
differential equations, we generalize the Huygens-Fresnel 
principle to obtain an integral formulation relating the 
values of the TV-photon wave function at a given time to 
its values on a fixed reference surface (usually a plane) 
at previous times. Familiar Fourier optics techniques can 
then be applied in order to engineer the propagation of 
entanglement and shape the photon state through the 
detection process. We illustrate this point using two ex- 
amples. Our approach generalizes previous works on the 
"ghost imaging" properties of photon pairs produced by 
parametric down conversion [3, 4] and is relevant to many 
applications of modern quantum optics, including quan- 
tum super-resolution imaging, quantum lithography, as 
well as spatial quantum communication. 

The very idea of using position wave functions to de- 
scribe the state of TV-photon systems relies on the ex- 
istence (but not uniqueness) of a photon position op- 
erator f , the Cartesian components of which are com- 
muting Hermitian operators satisfying [fk,Pi] = ifr$ki, 
p being the photon momentum [5, 6]. The eigenfunc- 
tions of f are transverse waves that can be interpreted as 
localized-photon states [7] . Any admissible single-photon 
wave function is obtained as a linear combination of these 
localized states. TV-photon wave functions are symmet- 
ric elements of the tensor product of TV single-particle 
Hilbert spaces. Because the definition of the position op- 
erator f is not unique, there is more than one way to 



assign a position wave function to a single photon. The 
most popular one is probably the so-called Bialynicki- 
Birula-Sipe wave function ip(r,t) = [ip + (r,t) ip-(r, t)] 
which has two vector components corresponding to pho- 
tons with positive and negative helicity [8, 9]. Each vec- 
tor component has a Fourier expansion that reads 



tf±(r,t)= / d 3 fc VhTce ± (k) /±(k) 



a i(k-r— kc t) 



(2tt) 3 / 2 



(1) 



where k = |k| and e±(k) are the unit circular polariza- 
tion vectors for photons propagating in the k-direction. 
Normalization is such that the complex coefficients f± 
satisfy Ylh=±J IA(k)| 2 = 1- This wave function 
transforms as an elementary object under Lorentz trans- 
formation and can be easily connected to Maxwell fields. 
In this Letter, we write the wave function is a slightly 
different (but equivalent) way that consists in summing 
both helicity components together: 



¥(r,t) =^+(r,f) + ^_(r,t). 



(2) 



This provides a vector representation instead of the bi- 
vector one [10]. Since ip+ and if>_ are orthogonally po- 
larized, they never mix: if ^ is given, and can be 
deduced. Therefore the information content in the vector 
function * is the same as in the bi- vector field ip. 

Working with wave functions ('first quantization" for- 
malism) is fully equivalent to working with quantum 
fields ("second quantization" formalism). Replacing 
the complex coefficients f± (k) by annihilation operators 
a±(k) in Eq. (1), the fundamental quantum field is found 
to be proportional to the positive frequency part of the 
electric field: 



h=± 



*(r,i) = d 3 k Vhk~ce h (k) a h (k) 



a i(k-r— kc t) 

(2tt) 3 / 2 



= -i V2ibE (+) (r,f). 



Note that 4? oc E( + ) holds information about both 
electric and magnetic field. Therefore, it provides a 
complete information about electromagnetic configura- 
tion. To show this explicitly, one decomposes * again 
into its helicity components ■0+ and and subtract 
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them: this yields B>( + ) = \J /io/2(ip + — i/>_), the pos- 
itive frequency part of the magnetic field. In the sec- 
ond quantization formalism, the state of a single-photon 
wave packet writes: |*) = J2h=± J d3fc A( k ) \ l Kh), 
where |lk,h) = a I(k)|0) and /±(k) are the same spec- 
tral amplitudes that appear in Eq. (1). The connec- 
tion between the first and second quantization formal- 
ism is given by the relation ^(r, t) = (0\4f(r, t)^} = 
-i y/2e^ (0|E(+)(r,t)|tf). Since ($|*(r,t)|tf) = for all 
|$) ? |0), we have ^*,(q')^i(q) = (*|*t,(g')^(g)|*) - 
2eo{E\7\q')E\ + \q)) for any pair of points q = (r,t) and 
q' = (r',t'), where the indexes € {x 7 y,z} 2 repre- 

sent Cartesian components. This relates the Bialynicki- 
Birula-Sipe wave-function to the usual first-order cor- 
relation functions of coherence theory. In particular, 
|*(r, t)\ 2 is proportional to the probability to detect the 
photon energy at point r at time t. This gives to the 
Bialynicki-Birula-Sipe wave function the usual interpre- 
tation of a probability amplitude to find the photon en- 
ergy at some position. 

The generalization to iV-photon states 



|*)= ]T J d 3 fcx... J d 3 k N 



h\,...,h 

//u,...,h JV (ki, . . • , kjv)|lki,/ni • • • j lk N ,/i N ) 

is straightforward. The connection between wave func- 
tions and fields is given by 

* il ... iN (q 1 ,...,q N ) = (-i) N (2e ) N ' 2 

x<0|£&W--4 +) (9i)l*> ( 3 ) 

and 

*3 , i...i' w (</i,"-,«5v)*ii...i w (9i,---,«/Jv) = i*o) N 
x (Elr\q[)...Elr\q , N )E^\q N )...E^(q 1 )). (4) 

Eq. (4) shows that any field correlation function of a 
iV-photon system can be computed as a product of two 
tensor elements of the A^-photon wave function. The best 
way to compute the propagation of the iV-photon wave 
function may depend on the situation. Often, the pho- 
tonic state is prepared in such a way that the wave func- 
tion is known at all times, but only on a specific surface 
E. Therefore, a diffraction theory of A^-photon states is 
needed. An important example is the generation of en- 
tangled photon pairs in a nonlinear crystal. In that case, 
E is the output face of the crystal and diffraction from 
that plane leads to phenomena the phenomena of "ghost 
imaging" [11]. 

To understand how A^-photon detection correlations 
spread in space and time, we first consider free-space 
propagation. We also make the simplifying assump- 
tion that we deal with paraxial states of light, in which 



case polarization does not change much during propaga- 
tion. Therefore, we drop the polarization-related indexes. 
Considering photons propagating along the z-axis, we 
use the Huygens-Fresnel principle [12] to express 
at some point as a function of its values on the reference 
plane E at z-coordinate (, 



\r-p\ 



(5) 



and inject this in Eq. (3). This yields 



*(ri,ti, . . .,r N ,t N ) 



(2ttc) 



dt 



-■■■^(pi,h- l -^,..., PN ,t N - 



|rjy-pjy| 



|ri - pi| ■ ■ ■ I r^v - Pn\ 

(6) 

We call this integral the generalized Huygens- 
Fresnel (GHF) principle for A^-photon wave func- 
tions. In Eqs. (5) and (6), Pj = (£j,VjX) 
(j e {1, ...,JV}) are points in the C-plane and 
pj- = (£,j,Tjj) arc their transverse components. In 
the optical domain, photons can usually be con- 
sidered as quasi-monochromatic. Therefore, the 
wave function can be written ^(ri, t\, . . . , rjv, tjv) = 

a(ri,ii,...,rjv,ijv)exp(-i27rc(|^ H h f^-)), where 

a(ri,ii, . . . , rjv,tjv) is a slowly varying function of time 
and Xj (j E {l,...,iV}) are the central wavelengths 
of the photons. Note that nothing prevents photons 
from having the same central wavelength or even being 
indistinguishable. Inserting that anzats in Eq. (6) and 
taking into account that a(ri, t\, . . . , rjy, tjy) is slowly 
varying in time, one obtains 



a(ri,ii, . . .,r N ,t N ) 



(-i) 



JV 



Ai...A 



N 



//dV^// d X, 



/ , l r i-Pi| , \rN-PN\, 
a (Pi,ti 7 • • • j Pn, in ) 

c c 

exp(iff|ri-pi|) exp(ij^\r N - pn^ 

|ri-pi| "' \r N - p N \ ' ^ 

Eq. (7) is only valid in free space. If propagation from pi 
to Yi is through an optical system, the free space propa- 
gator 



h fs {Yi,pi 



-i exp^x-l 1- *-^) 



A, 



\r% - pi 



(8) 



must be replaced by the appropriate one hifa, p{). With 
this generalization, Eq. (7) becomes 

a(r 1 ,t 1 ,...,r N ,t N ) = JJ d 2 p{ . . . j j d 2 p N 

i . l( r i,Pi) , l(r N ,PN), 
a (Pi,ti 7 • • • 7 Pn, tjv ) 

c c 

h\{r 1 ,p{)...h N (T N ,p N ), (9) 



3 



where Z(rj,p,) is the optical path length from pi to r^. 
Formula (9) assumes that there is only one optical path 
from pi to Vi. However, interferometers with arms hav- 
ing different path lengths can be placed between pi and 
rj. To take this into account, we generalize (9) in the 
following way: 

a(r 1 ,t 1 ,...,r N ,t N ) = J J d 2 pi...JJ d 2 p^ 

\- , , Mj^Pi) h N (r N ,p N ) 
a(p!,ti ,...,p N ,t N ) 

fei,...,fejv 



h[ kl) (r 1 ,p 1 )...h^ N '(r N ,p N ). (fO) 

The indexes h label the different paths from pi to r-j. 
We illustrate formula (10) using two examples. Con- 
sider the system in Fig. 1. The scheme is similar to the 
original "ghost imaging" scheme [11], with the difference 
that two diffraction masks (called Mi and M 2 ) are placed 
in the arms of a balanced Mach-Zcnhder interferometer. 
Co-propagating energy-degenerated (Ai = A 2 = A) en- 
tangled photon pairs are generated by parametric down- 
conversion in a nonlinear crystal and split using polariz- 
ing beam splitter (PBS). The generation process entan- 
gles the photons in position and time. The detector Db 
clicks whenever the photon going upward passes through 
Mi or M 2 and is detected on axis. The diffraction pattern 
of the second photon is detected by moving the point-like 
detector D a . Assuming a Gaussian pump and broadband 
phase-matching, the wave function amplitude at the out- 
put of the nonlinear crystal is well approximated by 



a(pi,h, p 2 ,t 2 ) 



exp(-^) exp(-%±^) 

(27TT 2 ) 1 /4 (2^ 2 )V2 

xS(pi-p^) S(t 1 -t 2 ), 



where T and S are the time and space rms-half-width of 
the pump pulse. The propagation in the optical system 
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FIG. 1: Interference effect in quantum imaging with two sep- 
arated diffraction masks. 



can be computed using formula (10), with 
-i exp [i^K - Pi I] 



hi(r a ,pi) = — 



|r« - Pi | 



V2\ 2 J J |r-p 2 | 
x cxp[i^|r b -r|] dV 
|r b — r| 

where M fc (r ± ) (fc e {f,2}) are the mask transfer func- 
tions. The GHF integral (10) simplifies if AS 2 > A(s + 
min(si,s 2 )) (i.e. the wavefront curvature of the gener- 
ated photons can be neglected) and if the far-field condi- 
tions d 2 /X <C S3 and d 2 /A <C (sq + s 2 ) apply, where d is 
the relevant length of the mask profile functions. Drop- 
ping constant and quadratic phase factors, one finds 



a(r a ,t a ,r b ,t b ) = exp 
2 



h - (so + S2)/c) 



exp [iQ]^M fe 



fc=i 



2ir ( x a 



AT 2 



S3 



5(t a - tb-r) 





( Va 








s 3/_ 



where L = (2s +si+.s 2 ), r = (s 1 -s 2 -s 3 )/c, i-; = (x l ,y i ) 
and Mfcfk- 1 ) = / M k (r^) expf-ik 1 - • r^] dr 1 - is the 
Fourier transform of the mask function Mfe(r _L ). The 
exponential factor indicates that the detection of a pho- 
ton in D u can only happen a propagation time (so + s 2 )/c 
after the photon has been created and within the initial 
time uncertainty T. The <5-factor shows that the differ- 
ence t in the detection times of D a and Db is just due 
to different propagation distances from the crystal to the 
detectors, as expected for time-entangled photons. The 
wave- front curvature of both photons are related through 
the quadratic phase Q = (n/X)[(x 2 + y 2 )/(2so + si + s 2 ) + 
( x b + Vb)l s z\- ^ detector Db, placed on axis (rjj- = 0), 
detects a photon a time % = t* , the wave function ampli- 
tude of the photon travelling towards D a is automatically 
projected on 



a(r a , t a ) = 6(t a — [t* + t]) exp (iQ a ) M k (2tt X 



fe=i 



al'HS;)' 



with Q a = (n/X)(x 2 a + y 2 )/(2s + s 1 + s 2 ). The system 
of Fig. 1 can be used to produce heralded single pho- 
tons with on-demand or adaptive spatial profile. One 
can mimic a complex modulation (amplitude and phase) 
using phase modulators at Mi and M 2 . Such a system 
would be useful in the context of earth-satellite quan- 
tum cryptography [13] if true single photons had to be 
used. Even with simple masks, such as the double slits 
in the inset of Fig. 1, non trivial shaping can be done. If 
2A = (5, changing the interferomctric phase <fi from to 7r 
doubles the spatial frequency of the interference fringes 
seen by D a . Single photons with very different spatial 
profiles can be created by tuning only one parameter. 
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As a second example, consider the system in Fig. 2. A 
nonlinear device produces photon triplets with time and 
space entanglement. (Such a device do not exist yet, but 
time-entangled triplets have been demonstrated recently 
[14]. A bulk version of that experiment would exhibit 
space entanglement as well.). We assume that, in the 
output plane of the device, the three-photon amplitude 



is 



a(pi,t 1 ,p 2l t 2 ,p 3 ,t 3 ) = 
x 5{pi - Ps) HP2 



cxp 



JL 

4T 2 



exp 



45 2 



(27rT 2 ) 1 /4 (2ns 2 ) 1 / 2 
pi) 6(h-t 3 ) 6(t 2 -t 3 ). 



We also assume that two photons have a common wave- 
length Ai, while the third one has a different one A 2 . 
They are separated using the dichroic beam splitter 
(DBS). In the Ai output of the DBS, we place a thin 
lens (focal /i) that images the output of the photon 
source in the plane of detector D\ with a magnification 
M = s 2 /(s + s\). The detector Db is placed on the 
optical axis. Using formula (9), one can calculate the 
three-photon amplitude in the detector planes. If a pho- 
ton is detected by Db at time f » , the wave function of the 
remaining Ai photons is projected on 

ex P I ~ iS'^M 1 

a{v a ,t a ,v' a ,t' a ) = 5{t a -{U+T)) S(t a -t' a )- ' 



(2T7S 2 ) 1 / 2 



5(x a x' a ) S(y a - y' a ) exp + I /M) 



h 2 (0,-r a /M). 



where r = (s 4 + s 2 — s 3 — s 4 )/c and h 2 (rb,p 3 ) is the 
propagator from the nonlinear device to the detector Db- 
In the plane of detector D a , the photons exhibit spatial 
bunching. The field is a linear superposition of local- 
ized two- photon Fock states of light. Due to the mag- 
nification factor M, the spatial extension of that linear 
superposition can be much larger than S. Such a quan- 
tum state of light is interesting in the context of coher- 
ent super-resolution imaging [15]: The object to be im- 
aged would be placed in the plane of detector D a in or- 
der to be illuminated with that special quantum state. 




3-photon 
device s 3 



Resolution enhancement only matters in optical systems 
in which geometrical aberration have been eliminated. 
Wavefront curvature must also be under control to avoid 
distortions due to non-isoplanatism [16]. The scheme of 
Fig. 2 makes it possible to control the wavefront cur- 
vature of the Ai photons by tailoring the detection of 
the A 2 photon. For instance, one can make the wave- 
front of Ai photons flat in the plane of D a by placing 
a lens (focal f 2 ) in the path to D 2 (see Fig. 2) and 
choosing S3, s 4 and f 2 such that 2(\ 2 / \\)M(M + 1) = 
s 2/((l/./2 — l/s4) _1 — (so + s 3))- This solution exists if 
the inequalities s 4 (s + s 3 )/(s + S3 + s 4 ) < f 2 < s 4 are 
satisfied. 

In summary, we developed a formalism that allows us 
to analyse many interesting issues related to the propa- 
gation of arbitrary number states of light using the wave 
function formalism. We derived a generalization of the 
Huygens-Fresnel principle that accounts for the propa- 
gation of field correlations (including entanglement) in 
space and time and showed how to applied it in practice. 
The formalism is very helpful to design sources of her- 
alded A-photons states with engineered spatial profile. 

This research was supported by the Interuniversity At- 
traction Poles program, Belgium Science Policy, under 
grant P6-10 and the Fonds de la Recherche Scientifique 
- FNRS (F.R.S.-FNRS, Belgium). 



FIG. 2: Generation of heralded linear superposition of local- 
ized two-photon states of light. 
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